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We study the magnetization dynamics in a nanocontact magnetic vortex oscillators as function of temperature. Low 
temperature experiments reveal that the dynamics at low and high currents differ qualitatively. At low currents, we 
excite a temperature independent standard oscillation mode, consisting of a gyrotropic motion of a vortex about the 
nanocontact in the free layer. Above a critical current, a sudden jump in the frequency is observed, which occurs with 
a substantial increase of the frequency versus current slope factor. Using micromagnetic simulation and analytical 
modeling, we associate this new regime to the creation of a vortex-antivortex pair in the pinned layer of the spin valve. 
This pair gives an additional perpendicular spin torque component that alters the free layer vortex dynamics, which can 
be quantitatively accounted for by an analytical model. 


I. INTRODUCTION 

Magnetic vortices are fundamental topological states 
in restricted geometries such as thin submicron dots or 
nanopillars. 1 2 For a certain range of aspect ratios, the micro- 
magnetic ground state is a vortex structure because the cir¬ 
cular configuration of the spins minimizes the stray dipolar 
fields. Because the norm of the local magnetization vector is 
conserved in strong ferromagnets due to a large exchange in¬ 
teraction, the magnetization at the core of the vortex tilts out 
of the film plane. 

Vortices can also be nucleated in extended ferromagnetic 
thin films in the nanocontact (NC) geometry 3 4 , in which it has 
been demonstrated that vortex manipulation over micrometer- 
scale distances is possible 5 . Spin torque effects appear when 
large current densities 10 12 A/m 2 ) are applied through the 
NC, which also results in significant Oersted-Ampere fields, 
e.g., 300 mT for 50 mA in 100 nm diameter NCs. Indeed the 
perpendicular component of this current flow (relative to the 
film plane) leads to a circulating Oersted-Ampere field akin 
to that generated by a cylindrical conductor. As such, a vor¬ 
tex state can appear by minimizing the Zeeman energy asso¬ 
ciated with the Oersted fields, in contrast to the case of con¬ 
fined geometries in which stray fields are minimized. Because 
the Zeeman interaction is proportional to the current, vortices 
only appear above a certain threshold or nucleation current 

C^nucl/ ■• 

However, processes involving vortex nucleation are sub¬ 
ject to conservation laws involving topological charged. The 
topology involved is described by the Skyrmion number, q = 
rjp/2, where rj is the vorticity which describes the curling 
magnetization of a vortex by p = +1 and p = -1 for an an¬ 
tivortex. p is the core polarity, which describes the orientation 
of the magnetization at the vortex core. Since a uniform state 
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has a total Skyrmion number (topological charge) of q = 0, 
the nucleation of a vortex ( q = p/2) must be accompanied 
by the nucleation of an antivortex with the same core polarity 
such that the total q remains zero 7 . 

The stability of a vortex-anti vortex (V-AV) pair in thin films 
strongly depends on the boundary conditions, i.e., on the mi¬ 
cromagnetic state at the system edges. It has been showrP 
that pair nucleation in the magnetic free layer in zero field 
is followed by the antivortex being expelled by the Oersted- 
Ampere field, resulting in steady-state oscillations of the vor¬ 
tex around the NCP. The resulting micromagnetic state of the 
free layer therefore resembles the well-known Landau state. 
In contrast, the presence of an in-plane magnetic field would 
favor an uniform magnetic state far from the NC. In this case, 
the antivortex would be bound to the vortex such that the uni¬ 
form state is preserved in the bulk of the film. This should 
equally be true for a ferromagnet exchange-biased by an anti- 
ferromagnet, where the internal field acting on the ferromag¬ 
net due to the exchange coupling also favors a uniform state. 

Previous studies® 1 have suggested that the two ferromag¬ 
netic layers of a nanocontacted spin-valve structure may con¬ 
tain a vortex state under certain conditions of the applied 
magnetic field and injected current. To date, a successful 
model to explain the new observed dynamics has been lack¬ 
ing. In this article, we address this question from experi¬ 
mental and theoretical perspectives. We present experimen¬ 
tal evidence of vortex-antivortex pair nucleation in the pinned 
layer of magnetic nanocontacts. At low bias currents, we ob¬ 
serve the usual free-layer vortex oscillations expected of such 
structures 12 . However, above a certain critical current, we 
detect the presence of a vortex-antivortex pair in the pinned 
layer through changes in the power spectra associated with the 
free-layer vortex oscillations. This critical current is strongly 
temperature-dependent, which is suggestive of a thermally- 
activated process. The experimental results are consistent 
with predictions based on rigid-vortex model. 
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FIG. 1. (Color online) Vortex oscillations in the free layer of the mag¬ 
netic nanocontact device, with uniform magnetization in the pinned 
CoFe layer. The vortex spirals out until a stationary orbit in a poten¬ 
tial provided by the Oersted-Ampere field. The vortex orbital motion 
leads to a time-varying voltage (top inset). Arrows represent the pro¬ 
jection of the average free layer magnetization underneath the NC 
area. 


II. EXPERIMENTAL RESULTS 

Figure 1 shows the experimental system studied. It con¬ 
sists of a metallic nanocontact fabricated on top of a spin- 
valve (SV) stack of width L = 17/im. The composition of 
the SV is IrMn(6)/Co 9 oFeio(4.5)/Cu(3.5)/Ni8oFe 2 o(5)/Pt(3), 
where the number in parentheses denote the layer thickness in 
nanometers. Details of the fabrication process are given else¬ 
where 13 . For the system studied here, the NC radius r nc is 75 
imP. To determine the free layer magnetization (Ni 80 Fe 2 o) 
and the Gilbert damping, ferromagnetic resonance experi¬ 
ments were performed from which we determined /xqM s = 1.1 
T and a=0.013, respectively. The electrical properties of 
the device were characterized by performing static magneto¬ 
transport measurements, which give a device resistance of 8.7 
D and a magnetoresistance of 25 mQ at room temperature. 
The CogoFeio layer acts as pinned layer since it is exchanged 
bias by an antiferromagnet (IrMn). Differential magnetoresis¬ 
tance curves have been also measured with a lock-in technique 
with a 10 fiA ac current and zero dc applied current for differ¬ 
ent applied temperature. This technique allow us to find out 
how the exchange bias field (//bias) varies from 80 K up to 420 
K, as it is shown in figure [3] 

Prior to electrical characterization of the nanocontact de¬ 
vice, we applied a magnetic field of « 115 mT along the easy 
axis in order to saturate the free layer magnetization and there¬ 
fore avoid to have any domain wall or vortex structure in the 
initial state. The device was measured at different temper¬ 
atures from 6 K to 300 K in a cryostat probe station under 
zero applied field. To characterize the magnetization dynam¬ 
ics, the corresponding high-frequency fluctuations in the giant 
magnetoresistance signal were measured with a spectrum an- 
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FIG. 2. Voltage power spectral densities (PSD) as function of cur¬ 
rent for different temperatures: (a) 6 K, (b) 40 K, (c) 160 K, and (d) 
200 K. 


alyzer after amplification. The dc current (I^ c ) applied to the 
nanocontact was ramped from 0 to 40 mA (upward scan) then 
back to 0 (downward scan), with electrons always flowing 
from the free to the pinned layer. The upward scans are used to 
determine the nucleation current I nuc \ of the free layer vortex, 
at which the voltage spectra exhibit a series of well-defined 
peaks representing the vortex gyration around the nanocon- 
tacP. 

Figure [2] shows the voltage power spectral density (PSD) 
as a function of the applied current, associated with oscilla¬ 
tions of the vortex for four different temperatures, 6, 40, 160, 
and 200 K. The spectra are measured from the highest cur¬ 
rent value (^ 40 mA) to the annihilation current (« 9 mA) of 
the vortex. Except for currents close to the annihilation of the 
vortex, where the dynamics is not well understood, a quasi- 
linear dependence of the oscillating frequency on current is 
observed. This result is consistent with a confining potential 
for the vortex dynamics that is determined by the Zeeman en¬ 
ergy associated with the Oersted-Ampere field, in line with 
previous observations on other spin-valve compositions®^ 
This ’’standard” oscillation mode corresponds to a vortex that 
orbits around the nanocontact in the free layer (Fig. 1). 

However, a different dynamical behavior is seen for cur¬ 
rents above I& ~ 30 mA at different temperatures. At this 
critical current (I cr it), and in the range of a few mA above this 
value, power spectra of the vortex oscillations exhibit a bi- 
modal character. We interpret this as a signature of thermally- 
driven hopping between two distinct dynamical modes with 
different frequencies, which we label hereafter as the upper 
mode (UM) and lower mode (LM). Above this critical cur¬ 
rent, the LM frequency branch represents a continuation of 
the standard mode, while the UM branch occurs at a higher 
frequency with a different slope compared with the LM mode. 
While we have not obtained time-resolved measurements that 
confirm the hopping between the UM and LM, a signature of 
the hopping is present in the PSD shown in Figs 2(a)-(d). If 
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FIG. 3. Critical current, nucleation current, and exchange bias field 
as a function of temperature. 


nanocontact region. This has two consequences on the vortex 
dynamics. First, it leads to a coupling between the gyrating 
vortex in the free layer and the vortex-antivortex pair through 
the dipolar interaction, which could lead to an additional term 
in the confining potential for the vortex. Second, and more 
importantly, the core magnetization of the vortex-antivortex 
pair leads to a perpendicular-to-plane component for the spin- 
polarized current flowing between the free and pinned lay¬ 
ers, which can generate additional spin torques for the vortex 
dynamics. The latter may explain the frequency jump at the 
threshold / cr i t along with the different slope df /dl observed 
for the UM. In this light, the existence of two oscillation 
modes suggests that the vortex-antivortex pair in the pinned 
layer has a finite lifetime, where the hopping is due to the 
repeated nucleation and annihilation of the vortex-antivortex 
pair due to thermal fluctuations. 


a vortex is hopping between two oscillation modes, UM and 
LM, it spends some time tuM and £lm in each mode. A clear 
feature of hopping is a signal that should appear at frequen¬ 
cies related with the inverse time that the vortex spends to go 
from one mode to the other one. The observed frequency at 
which this mode appears is « 50 MHz. At high current there 
appears a low frequency shoulder that only survives as long 
as the UM and LM modes are present. We have gathered in 
Table I the frequency jumps /up-/lm at 7 cr i t and their slopes 
(df /dl) at different experimental temperatures. The depen¬ 
dence of J nuc i and J cr i t as function of the temperature is shown 
in Fig. m it is worth noticing these two characteristic currents 
follow the same trend, which suggests they involve similar 
physical processes. While / nuc i represents the threshold cur¬ 
rent to nucleate a vortex state in the free layer, the value J cr i t 
represents a threshold current associated with some dynamics 
in the pinned layer (PL). The observed decrease in power (not 
shown) when the current is increased is in agreement with this 
assertion. 

The exchange bias field acting on the PL decreases as 
the temperature is increased. Note that the cryostat temper¬ 
ature in our experiment does not take into account the Joule 
heating resulting from the current flow through the nanocon¬ 
tact device. Recently, it has been reported how the current 
density is distributed along the SV stack after passing trough 
the NCP. This study allowed us to simulate the temperature 
profile across the SV thickness underneath the NCP. For cur¬ 
rents of about 50 mA, we found the temperature increase to be 
between 150 K and 200 K in the vicinity of the NC. For exam¬ 
ple, at room temperature the H bias is about 30 mT whereas for 
an injected current of 50 mA it would decreases down to 22 
mT. Therefore, one consequence of Joule heating is to reduce 
the bias field acting on the PL in the nanocontact region. In 
line with this reasoning, nucleation of some magnetic struc¬ 
ture with an out-of-plane magnetization may occur in the PL 
underneath the NC. This is discussed in more detail in Section 
IV. 

The presence of a V-AV pair in the pinned layer would 
give rise to an out-of-plane component of magnetization in the 


III. MODEL FOR VORTEX-ANTIVORTEX 
NUCLEATION 

In order to quantify the scenario involving nucleation and 
annihilation of the vortex-antivortex pair described in the pre¬ 
vious section, we have extended the Thiele formalism for de¬ 
scribing the nanocontact vortex oscillations by accounting for 
the presence of the vortex-antivortex pair. We first examine 
the nucleation problem without and with spin transfer torques, 
and then provide a description of the vortex dynamics in the 
presence of the vortex-antivortex pair. 


A. Vortex-antivortex nucleation in the pinned layer in 
the absence of STT 

As shown previousl)®, the onset of vortex oscillations in 
the free layer takes place after the following sequence: (1) 
the initial uniformly-magnetized state is distorted by the large 
Oersted-Ampere field; (2) nucleation of a vortex-antivortex 
pair occurs after this distortion becomes irreversible; (3) ex¬ 
pulsion of the antivortex away from the nanocontact region. 
A similar process is expected for vortex-antivortex pair nucle¬ 
ation in the pinned layer, but the key difference is the presence 
of the exchange bias field i?bias that will limit the separation 
distance between the vortex and the antivortex. The equilib¬ 
rium separation will therefore be determined by balancing the 
competing attractive and repulsive forces. 

To describe the relevant energies associated with nucle¬ 
ation, we use the rigid vortex model to describe the vortex- 
antivortex pair. This formalism allows us to express the rel¬ 
evant energies in terms of the positions of the vortex and an¬ 
tivortex cores by using a suitable ansatz for the spin structure 
at the cores. To simplify the integrals for the energies, we as¬ 
sume that the vortex is centered about the nanocontact, X v = 
(0,0), while the antivortex is situated at X=(X, Y). The spa¬ 
tial distribution of the core magnetization in polar coordinates 
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(0 and <J>) is taken to be, 


$(x,y;X,Y) = rj tan 1 ( 1 ) 

where 77 = ±1 is the vorticity. By using this ansatz , the total 
magnetic energy of the pinned layer, 

E = E e xc H - Eoe H - Eebi (2) 

which represents the exchange interaction, Oersted field Zee- 
man energy, and the exchange bias interaction, respectively. 
The expressions for the different energy terms are as follows, 

E e xc = A p id p i j d 2 r (V $) 2 (3) 

E 0e = -HoHi • M'f j dV cos ($ 0e - $ c ) (4) 


E eb = -fj, 0 H bias ■ Mf ■ d p i J d 2 r cos ($ s ) (5) 

which can be simplified to give 

Eexc — “InApidpi In (6) 

r c 

E 0e = -Mo Mf l H r (nd p ir nc Xy) In — , (7) 

J* nc _ 

Eeb = MO HbiaeMfd? 1 *. (8) 


Here, A v \d v \ is the exchange stiffness (0.4 eV), H\ = 
\I\ /47rr nc represents the Oersted-Ampere field, <J>oe and 
describe the magnetic texture of the Oersted-Ampere field and 
the antivortex respectively, /i 0 M s pl =1.56 T, d v \ represents the 
thickness of the pinned layer (see Fig. 1), and k represents the 
effective surface area covered by the vortex-antivortex pair. 
For a vortex-antivortex separation in the range of 0 to tens of 
nm, this surface can be approximated by a disk joining the two 
cores and can be expressed as \\Xy || 2 . 

Figure |4ja) illustrates how the barrier energy related to pair 
nucleation depends on the V-AV separation when a dc current 
of 32 mA (I cr it) is applied. When the V-AV pair is present, the 
energy minimum corresponds to a given separation between 
the cores, which depends on the magnitude of the exchange 
bias field. Note that the Oersted-Ampere field was calculated 
considering the NC as an semi-infinite cylinder. This is known 
15 to slightly overestimate the Oersted-Ampere field, hence 
the corresponding energy term (eq. [7]); with our material pa¬ 
rameter, we will see that this can lead to an underestimation 
of the separation between the cores of the vortex and the an¬ 
tivortex. 

The model predicts that different magnetic states of the PL 
can appear depending on the magnitude of the field H bi as . For 


Ffbias > 40 mT, the energy as function of the V-AV separa¬ 
tion distance would exhibit a minimum in the range of the 
exchange length, 7 exc « 6 nm. Because this value is much 
smaller than the spatial extension of the vortex-antivortex pair, 
it suggests that most likely ground state in a pinned layer 
with T/bias > 40 mT would be simply a distortion of its 
magnetization underneath the NC. In the other limit (’’un¬ 
pinned” layer with JTbias — 0 mT), the energy minimum 
would correspond to a V-AV separation distance that falls out 
of the range showed in Figure |4ja). The preferred state in 
this case would thus involve the expulsion of the antivortex 
outside of the nanocontact region with the vortex remaining 
centered about the nanocontact. For the intermediate cases, 
0 < Tfbias < 40 mT, the exchange bias field leads to an en¬ 
ergy minimum at a well-defined separation distance between 
the vortex-antivortex pair; this distance has a strong depen¬ 
dence on the magnitude of the exchange bias field. We note 
that under the highest value of the applied current we used 
in our study (40 mA), the Joule heating leads to a tempera¬ 
ture increase of 120 kP above room temperature, which re¬ 
sults in a bias field of 7/bias ~ 22 mT; JTbias vanishes at about 
600 K in our samples. We conjecture that the separation dis¬ 
tance of the vortex-antivortex core remains sufficiently small 
in our experiment such that spin torques, thermal effects, or 
both combined can lead to thermally-activated pair annihila¬ 
tion. This provides a mechanism for the intermittence of the 
LM and UM we observe in the power spectra. 

To shed further light on this scenario, we performed micro¬ 
magnetics simulations on the nucleation process in the pinned 
layer using the Mum AX code 17 . The simulations were per¬ 
formed at zero temperature without the spin-torque terms due 
to the currents flowing perpendicular to the film plane. The 
material parameters used were M s = 1.260 kA/m, A = 19 
pj/m, and a = 0.013. The simulated region was a rectangular 
volume with dimensions of 1280 nm x 1280 nm x 5 nm that 
was discretized using 512 x 512 x 1 finite-difference cells. 
The spatial distribution of the Oersted-Ampere field was com¬ 
puted with full 3D finite-element simulations (COMSOLjEl, 
which was then included into the micromagnetics simulations. 

Results from the micromagnetics simulations are shown in 
Fig.[4jb)-(e). Our simulations confirm that a stable separation 
distance between the vortex-antivortex pair is possible for a 
given value of the applied current and exchange bias field. We 
notice that the main difference between the analytical model 
and the simulations is the fact that the antivortex is found to 
be stabilized outside the nanocontact area while in the ana¬ 
lytical model the antivortex remains inside the nanocontact. 
While there are quantitative differences in the separation dis¬ 
tance, the simulation results supports the physical picture that 
underlies the analytical description. 


B. Vortex-anti vortex nucleation in the pinned layer in 
the presence of STT 

The previous simulations neglected spin transfer torques 
(STT) in the nucleation process of the V-AV pair in the pinned 
layer. In the present step, we now take into account the STT 
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FIG. 4. (Color online) (a) Topological configuration of the PL for Id C =32 mA. The dashed lines separate the regions without pair, with a 
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arising from the intralayer spin currents with a spin polariza¬ 
tion constant P = 0.2; the other system parameters remain 
unchanged and are still meant to describe a single isolated 
pinned layer. Results from these simulations are given in Fig. 
5. As expected, the addition of the STT results in the vortex 
spiraling out of the edge of the NC. The simulations suggest 
that the vortex never reaches a stationary orbit because it anni¬ 
hilates with the antivortex first, resulting in a distorted micro- 
magnetic state in the pinned layer. After this annihilation, the 
process of nucleation can repeat itself. This scenario accounts 
for the intermittent hopping between the two vortex oscilla¬ 
tion modes (see Fig. [2]). The UM is likely to correspond to the 
existence of the pair while the LM is likely to correspond to a 
distortion of the pinned layer (without any out-of-plane com¬ 
ponent of the magnetization). If the current is considerably 
larger than J cr i t , e.g., 50 mA, the antivortex is expelled from 
the pinned layer and the vortex remaining in the pinned layer 
starts to perform full rotations around the NC along a station¬ 
ary orbit. At these very large currents, we would thus be in 
a situation in which the dynamics inside the pinned layer and 
inside the free layer are qualitatively similar. 


C. Influence of a V-AV pair in the PL onto the 
oscillating vortex of the free layer 

The analytical theory presented here is based on previous 
studies in which it was shown that the STT torques from 
the interlayer spin currents alone cannot drive self-sustained 
oscillations under zero field. Instead, it was shown that this 
is the in-plane (i.e. intralayer) current that drives the vortex 
motion. However this was for a static and perfectly in-plane 
magnetized pinned layer. 

Here the presence of the V-AV pair in the PL induces a 
slight out-of-plane tilt of the spin polarization of the CPP cur¬ 


rent. To describe the free layer dynamics, we need to take into 
account this slight out-of-plane tilt of the spin polarization of 
the CPP current P± = p± f dV sin 2 ©V^>, where (0, T>) rep¬ 
resent the magnetization orientation in polar coordinates. In 
our specific case, the prefactor of the integral, p± is calculated 
as the ratio between the vortex core radius and the radius of the 
NC giving a value for p_\_ « 0.02. The value of p± depends 
on the material parameters in the sense that the vortex core 
radius scales with the exchange length (Z exc ) of the magnetic 
material. Therefore, the larger the saturation magnetization 
M s , the lower the value of Z exc will be and therefore p± is ex¬ 
pected to decrease accordingly. The impact of the perpendic¬ 
ular component of spin transfer torque acting over the vortex 
in the free layer follows Thiele’s approach 20 . It includes non 
conservative torques and it is possible to describe the vortex 
motion around the NC. To describe the magnetization orien¬ 
tation we use polar coordinates 0(f) and <F(f). X = (X, Y ) 
represents the vortex position in the free layer. The equation 
of motion can be expressed as follows, 


Gx 


f 


+D 




+ «IP^= 0 . 

dX 


(9) 

where U z represents the Zeeman energy, a represents the spin- 
torque efficiency, and u(f) the spin-current drift velocity, a 
and (3 represent the damping constant and the nonadiabatic 
constant, respectively. By solving Eq. 6, using R 0 exp(i^) = 
XqPiYq, we find the following coupled differential equations, 


Gd t R + aDRBt -p + Gu = aIP± sin 2 e 0 , (10) 

jR 


GRd t p - aDd t RY (3Du = (11) 

where <4 and M s , represent the thickness and the saturation 
magnetization of the free layer, respectively. To determine the 
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FIG. 5. (Color online) Micromagnetic simulations showing the time evolution of the vortex-antivortex pair when a constant current Idc of 40 
mA is applied, for different exchange bias bias field, (a) /xo Tobias = 40 mT and (b) /xoi7bias = 18 mT. The intralayer STT is included in this 
simulation. 


radius of the stationary orbit of the vortex we set d t R = 0, 
which gives 


Ro = 


aIa 2 nc G 2 (1 + \p _l 

T^core ) f dU < 


aD 


dR 


( 12 ) 


This leads us to classify the FL dynamics into two regimes, 
depending on the absence (p±= 0) or presence (p± / 0) of 
a V-AV pair in the PL. The effect of the PL vortex is to in¬ 
crease the radius of the orbit of the FL vortex. Micromag¬ 


netic modeling, in agreement with Eq. 12 yields Ro=l 10 nm 
when there is no pair in the PL, and 110 + 10n core when cores 
are present in the PL underneath the NC. The corresponding 
oscillation frequencies of the FL vortex can be estimated if 
we assume that the Gilbert damping a and nonadiabatic spin 
torque parameters /3 have similar magnitudes 2122 . The fre¬ 
quency tunability, defined as the slope of the frequency versus 
current relation, then depends on the number of cores (n core ) 
underneath the NC, 


duiFL 


1 


dl GR 0 V dR 


dU z 


n c 


, aD 


G 


-crp_L. 


(13) 


The agreement between Eq. 13 and the experiments is illus¬ 
trated in Figure [6] where we show the dependence of the os¬ 
cillation frequency versus the applied current. Table I sum¬ 
marizes the theoretical and experimental values for the slope 
(doj osc /dl) and the frequency jump for different temperatures. 

The validity of the model can be checked by comparing the 
analytical and experimental slopes. By using the experimental 
values of /xoM s , d v \, a , and L, and by assuming a spin polar¬ 
ization of P = 0.5 and a radius of the vortex core r core =10 nm, 
we find a theoretical value for duj osc /dI =4 MHz/mA in the 
low frequency regime, which matches with the observed slope 
experimentally of « 3.6 MHz/mA from Eq. (10). The jump 
in frequency can be calculated theoretically and accounts for 
a jump of 36 MHz similar to the observed one of « 37 MHz. 
The larger slope in the high frequency regime (upper mode) 


TABLE I. Frequency jumps and slopes versus current at different 
temperatures. LM and UM stand for lower and upper modes. n cor e 
represents the number of vortex or antivortex cores in the pinned 
layer underneath the nanocontact. 


T 

(K) 

flJM — /LM at /crit 

(MHz) 

df/dl (LM) 
(MHz/mA) 

df/dl (UM) 
(MHz/mA) 

Theory (0) 

«36 

« 4 

~ 4+2 77-core 

6 

40 

4 

8.7 

80 

36 

2.8 

9 

120 

44 

3.4 

5.8 

160 

40 

3.6 

5.7 

200 

28 

4.3 
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results from the joint contributions of the Zeeman potential 
and perpendicular spin torque. This gives a larger slope which 
again fits well with the experimental one (Table I). 


IV. CONCLUSIONS 

In summary, we have shown that the nucleation of a vortex- 
antivortex pair in the pinned layer of a spin valve nanocon¬ 
tact oscillator can lead to distinct changes to the oscillatory 
dynamics of the vortex in the free layer. The pair leads to an 
additional spin torque term related to currents flowing perpen¬ 
dicular to the film plane, which results in a frequency jump 
along with a different frequency tunability for the free layer 
vortex gyration. The temperature dependence of this effect 
is related to variations in the exchange bias field acting on 
the pinned layer, which determines the equilibrium separation 
distance of the vortex-antivortex pair. Intermittent modes in 
the free layer dynamics are attributed to thermally-activated 
pair nucleation and annihilation in the pinned layer. Analyt¬ 
ical modeling and micromagnetics simulations are shown to 
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FIG. 6. (Color online) (a) Sketch of the standard vortex mode and 
its new dynamics when there is a vortex-antivortex pair in the pinned 
layer. The extended colorbar below the sketch represents the mode 
at which the vortex is oscillating depending on the applied dc cur¬ 
rent. Panels (b), (c), (d), (e) show the mode frequency at different 
temperatures, 80, 120, 160 and 200 K, respectively. 


give good agreement with experiments. 
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